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Abstract 



Loss processes that remove particles from an atom trap leave holes behind 
in the single particle distribution if the trapped gas is a degenerate fermion 
system. The appearance of holes increases the temperature and we show that 
the heating is (i) significant if the initial temperature is well below the Fermi 
temperature Tp, and (ii) increases the temperature to T > Tp/4 after half of 
the system's lifetime, regardless of the initial temperature. The hole heating 
has important consequences for the prospect of observing Cooper-pairing in 
atom traps. 
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Insofar as the ultra-cold atoms behave as an equilibrium system, the neutral atom trap 
provides a laboratory for the study of low-temperature phenomena. Recent experiments with 
atomic Bose-Einstein condensates (BECs) |1J further highlight this connection by reporting 
genuine superfluid behavior |2[]. Similarly, the observation of fermion pairing in the gas 
phase would be of equal, or even greater interest |§, and the quest for atom-trap fermion 
superfluidity has begun in earnest. In this Letter, we discuss the effects of particle- loss to 
which the metastable ultra-cold atom gases are prone. In addition to limiting the system's 
lifetime, t l , the loss processes leave behind holes in the single particle distribution of a 
degenerate fermion gas (of temperature T significantly lower than the Fermi-temperature 
Tp). The fermions heat in the subsequent relaxation and, without additional cooling, the 
system's temperature doubles from its initial value T within a time t<i ~ 20 [T /Tp] 2 r^. After 
t ~ 0.5t£ the fermion gas heats up to T > 0.25Tp, regardless of the initial temperature. 
The competition of superfluid formation with 'Fermi- hole heating' leads to stringent lower- 
bounds for the strength of the inter-atomic attraction that would pair the fermions. 

Introduction In atom traps, inelastic collisions impart kinetic energy to the scattering 
products, which subsequently leave the trap. Scattering with background atoms and inelas- 
tic two-and three-body collisions decrease the local particle density, n, as n = —•yn, n = 
—/3n 2 ,h = —an 3 , respectively. Although two and three-body scattering is suppressed in a 
single fermion gas, so are the inter-atomic interactions that relax the system to thermal equi- 
librium. By trapping 40 i^-atoms in two different internal states D. Jin's group circumvented 
the problem and succeeded in cooling fermion atoms evaporatively (i.e. by removing the 
highest energy atoms) to degeneracy Q. Most recently, R. Hulet's group at Rice University 
H and C. Salomon's group at the ENS, Paris 0, cooled fermion atoms ( a Li) by bringing 
them in thermal contact with ultra-cold ( 7 Li)-bosons. The lowest reported temperatures 
H are T ~ 0.257V, and the fermion mixture cooling efforts appear to have encountered a 
limit in this temperature range M. Calculations suggest that Pauli-blocking plays a role 
0] but this effect alone does not seem to explain the limit. An alternative motivation for 
creating fermion gas mixtures is the prospect of s-wave Cooper-pairing. This scheme - 
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the most likely mechanism to achieve fermion superfluidity at realistic atom trap densities, 
n ~ 10 12 — 10 15 cm~ 3 - requires the s-wave interaction between two types of distinguishable 
atoms to be attractive (i.e. described by a negative scattering length, a < 0). Whether 
the atom mixtures are created to equilibrate the system or to form Cooper-pairs, the unlike 
fermion interactions introduce loss processes. 

We consider a weakly interacting gas mixture of the same fermion isotope in two different 
internal states, 1 and 2. The optimal mixture contains equal particle densities n\ = n<i = n, 
corresponding to equal particle numbers, Ni = N 2 = N ||. The usual trapped neutral 
atom gases are dilute in the sense that \a\/r << 1 where r denotes the average inter-particle 
distance, r = n" 1 / 3 (typically, \a\ ~ lnm whereas r ~ 1. — 0.1/xm). In contrast, the quest 
for fermion superfluidity necessitates an unusually strong neutral atom attraction. This 
condition follows from the critical temperature T c for s-wave Cooper-pairing: in the described 
homogeneous fermion mixture of Fermi-momentum kp = (67r 2 ) 1 / 3 r _1 rs 3.9r _1 , Fermi-energy 
tp = h 2 k F /2m, and Fermi-temperature Tp = epjkp (where kp is the Boltzmann constant) 
T c equals 



T c « 0.3T F exp 2 M«i =0.3TVexp u ~ . (1) 

With Tp w S^fiKA-^n/no) 2 / 3 , where no is a realistic reference density, uq = 10 cm , 
and A the mass number, we find that the prefactor of (|l|) is accessible to present day atom 
trap technology, but with \a\/r ~ 10 -3 , T c w 0.3Tp exp(— 400). 

Lifetime restrictions Even if such exponentially low temperatures can be reached, the 
time scale on which the superfluid forms, r/ orm , exceeds any realistic lifetime Tp. This for- 
mation time, Tf orm , depends on the Fermi-time scale, Tp ermi = h/ep ~ 2.09 fxsecA(n/ n )~ 2 / 3 
and takes on the form flT0| 
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^■T Fer . mi exp(0.8r/|a|). (2) 

Z7T 



The requirement that the superfluid forms before the gas is depleted, Tf orm < tl, con- 
strains the scattering length according to \a\/r > 0.8/[Zn(27r/10) + ln(r L / T Fermi )] % 



0.35/logi (r L /r Fermi ) (assuming r L >> 10 2 T Fermi ). For typical values, TF ermi ~ 10 — 100/isec 
and t l ~ 100 seconds, logi (r L /r Fermi ) ~ 6 — 7, we find that the scattering length has to 
exceed 5% of the average inter-particle distance, \a\/r > 0.05 — 0.06. Atom gas Cooper- 
pairing thus requires an effective interaction that is unusually strong for neutral atoms. Such 
interaction can be realized by using atoms that form a virtual state in the binary scatter- 
ing, or by using external fields that alter the atom-atom interactions. On the other hand, 
when the magnitude of the negative scattering length exceeds \a\ > 0.48r (corresponding to 
kp\a\ > 37r/5, |§), the pressure turns negative and the gas collapses. This requirement of 
mechanical stability limits the critical Cooper-pairing temperature to T c < T c ^ max pa O.lTp. 

Hole heating Next, we illustrate the heating caused by the loss-induced creation of fermi- 
holes. We consider a homogeneous fermion mixture that is 'normal' (i.e. not superfluid) 
and initially at zero temperature. The single particle distributions are filled fermi-spheres 
in momentum space, corresponding to occupation numbers n^k — 1 (j — 1, 2) if k < k F , 
n j,k — if k > k F . The loss processes 'perforate' the fermi-spheres, creating holes that bring 
the system into a state that is not the ground-state (which has filled fermi-spheres of reduced 
radius). In due course, the system relaxes to its thermal equilibrium at finite temperature. 
The scattering processes that relax the system are interesting in their own right. A hole in 
Fermi-sphere 1 is filled by a particle of initially higher energy. This fermion 1 can change its 
momentum by interacting with a fermion 2 that is thereby promoted to an energy above the 
Fermi-level. Such particle-hole scattering, akin to traditional Auger scattering, continuously 
produces particles with energies up to 2e F , energies that are anomalously high from the 
perspective of a thermal distribution. The resulting 'high-energy' fermions are scattered to 
lower energy states (or 'cooled'). 

The importance of hole heating follows from the rate of temperature increase which 
depends on the specific heat per particle, cy = (n 2 /2)/c B [/c B T/eF] and on the rate with 
which particles are removed, 



where the system's lifetime, tl, can depend on the density (t^ 1 = 7, /3n, an 2 for background, 
two- and three-body loss respectively) The removal of a single particle of momentum 
k lowers the energy of iV fermions at temperature T, E(N,T), by its kinetic energy but 
increases the energy relative to that of the remaining particles, N — 1 in total, at the same 
temperature, 
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where we assumed T << Tp so that the chemical potential is approximately equal to the 
Fermi-energy: E(N,T) — E(N — 1,T) rs tp. Therefore, upon removal of a particle of 
momentum k, the system 'gains' [ep — h 2 k 2 /2m] in excess to the energy of the remaining 
particles at the same temperature. The probability to create a hole in the (J, k)-state after 



a short time interval At is equal to —At hj^ = — %,k- The total excess energy gathered 
by the system during At is the (j, k)- excess energy multiplied by the probability for hole 
creation, summed over all (j, k)-states: 
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The corresponding temperature increase AT = AE /[2Ncy] gives a rate, T = AT/ At, equal 
to 
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The solution to Eq.flf]), parametrized by the initial temperature T and the 'temperature 
doubling time', r 2 , 

157T 2 



T = Tjl + - ;r 2 
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shows that the system doubles its temperature in a time that is proportional to the square 
of the Fermi to initial temperature ratio, r 2 ~ 20[T /T f ] 2 t l . Thus, fermions brought to 1% 
of their Fermi-temperature double their temperature in ~ 0.2% of their lifetime. 



The above derivation assumes a constant loss-rate, whereas a change in the density alters 
this rate on a long enough time scale. Also, the use of the specific heat, a statistical quantity, 
may give the impression that the validity of the above derivation requires the system to 
remain in thermal equilibrium. To address both issues we determine the temperature from 
the observation that the energy-per-particle, e, of a homogeneous fermion gas or fermion 
gas mixture, is not affected by the particle loss or by the subsequent scattering that 
brings the system to thermal equilibrium. In the Fermi-degenerate regime, we know that 



e(N,T) « (3/5)e F (N) + (vr 2 /4)[(k 2 B T 2 )/e F (N)] } correct up to order [T/T F f. Equating 



e(N, T) to its initial value e(No, To), we solve for T. Denoting the initial Fermi-temperature 



by T F>0 , using that [e F (N )/e F (N)] = [iV /iV] 2 / 3 , and [T F /T F>0 ] = [N/N ] 2 / 3 , we obtain 



N(t)/No = exp(— t/rL), [1 + t/ri]^ 1 , [1 + 2£/tx]~ 1 / 2 for loss that is predominantly caused 
by background, two- and three-body scattering respectively. In Fig.(l), we graph the 
time-dependent temperature to Fermi-temperature ratio for T = O.OlTp. The short-time 
behavior, t < O.Itl, is indeed well described by Eq. ([/]). On the scale of the graph 
(t ~ r L , T ~ O.lTp) the plot is 'universal' in the sense that any temperature curve of 
lower initial value is indistinguishable from that shown (described by the To — > 0-limit, 
T — > T F fl\J (12/ hn 2 )\(N / Nq)~ 2 1 3 — 1]). After t ~ O.Itl, the temperature increase slows down 
and T reaches 25% — 30% of the Fermi-temperature by 'middle age', t ~ 0.5rj,. At later 
times, the fermion system becomes classical (invalidating the above expression of e(N,T)) 
and lim^oo T(t) = 0.4T Fj0 , which follows from e dassica i = (3/2)k B T ps (3/5)e F (iV ). The 
above result does not depend on thermal equilibrium, provided we interpret temperature as 
a measure of the difference of the system's energy with its ground state value. 

Superfluid formation and heating as competing effects Regardless of the initial temper- 






To obtain a time-dependent expression, we substitute the fraction of remaining particles by 
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ature, the fermions heat up to at least 10% of the Fermi-temperature after only 7% of the 
system's lifetime. The temperature T = O.lTp is the highest possible critical temperature 
for fermion superfluidity, T c ^ max ph O.lTp, therefore, the maximal time the fermion super- 
fluid has to form is 0.07tl, not tl- It may, in fact, be difficult to realize T C)Tnax and we 
consider the more general constraints on the interaction from the requirement of forming 
the Cooper-paired superfluid before the system heats up above the critical temperature: 

T(t = T form ) < T c . (9) 



Because the superfluid formation time cannot exceed 0.07r L , T yl + 3t/r 2 (|7|) adequately 
describes the temperature, and we find from (§) 

T2 > ( w - 1 ■ < 10 > 

Substituting r 2 and Tf orm in terms of r L (0) and T Fermi @, the inequality (|H]) yields 

(J\ (T F /T C ) 4 

i-(T /T c y Fermi - [ } 

If the experimentalist succeeds in bringing the initial temperature of the system significantly 
below the critical temperature (e.g. To < 0.3T C ), then the denominator in flTT| ) can be 
replaced by 1 and the condition becomes independent of T . Assuming T < 0.3T C and 
substituting T c /Tp by 0.3 exp(— 0.4r/|a|), we finally obtain 

\a\ 0.4 

> 



r ln[0.3/(4/57r 3 ) 1 /4] + (1/4) \n(r L / r Fermi ) 
iogio{T L /T Fermi ) 

where the last approximation is valid if tl > 10 4 r^ erm j. Note that the minimal value of the 
dimensionless interaction strength \a\/r is approximately twice as large as from the condition 
Tform > t l- Thus, if T L /r Fermi ~ 10 6 — 10 7 , the minimal value the negative scattering length 
has to take on to form Cooper-pairs is, \a\/r > 0.125 — 0.106. 

Equation (([J) tacitly assumes that Fermi-holes relax faster than the superfluid can form, 
an assumption we now justify by calculating the lifetime of a hole. Specifically, we consider 



a fermion 1 hole in the center of the Fermi-sphere (l,k = 0), and we determine its inverse 
lifetime t^ =0 in a Fermi-Golden rule calculation as the rate at which the hole 'fills up" by 
scattering processes (1, q) + (2, p) — > (1, k = 0) + (2, p + q): 

2 

^ni jq 72 2 ,p(l — ™2,p+q) <5 '(E 'initial ~ Ef ina i) (13) 

p>q 

where Q denotes the volume in box normalization. With the (1, k = 0) final state, the energy 
density reduces to S(Ei n i t i a i — Ef ina i) — > (m / 'h 2 pq) S (cos 9) , where 9 denotes the angle between 
q and p. The (q, p)-integral with zero-temperature Fermi Dirac occupation numbers yields 
r Mc=o = (4vra 2 )n ) • Remarkably, t^I_ =0 equals the classical estimate for the collision 
rate of a zero velocity particle in the presence of the distinguishable particles, where 
Tcoii = [(crnv)]^ 1 = .207 TFermi(r / \ a\ ) 2 . Replacing the Pauli-blocking factor (1 — ^2,p+q) in 
(|T3|) by 1 reduces the lifetime to half its value. Therefore, whereas Pauli-blocking increases 
the lifetime of holes and particles drastically near the Fermi-surface, its effect on the hole 
lifetime is moderate near the center of the Fermi-spehere where r^k ~ t co u, which is also 
the time scale on which the high energy (2, p + q)-fermions scatter. Finally, Tf orm /r co u « 
7.7(|a|/r) 2 exp(0.8^|), so that at its minimum (|a| = 0.4r) T form /T coU w 9.5 and T form /T coU 
quickly increases as \a\ is reduced. The implication is that that there is a region inside 
the Fermi-sphere where holes are short-lived on the time scale of superfluid formation. The 
decay of those holes produces 'high energy' fermions that start heating the system while the 
superfluid forms. 

While the treatment (|9|)-(|i"2"|) certainly simplifies much of the intricate dynamics, we can 
trust the accuracy of the final constraint ( |12"D to at least 15%. Consider, for instance, the 
possibility that the non-equilibrium nature of the formation process, or the change of cy as 
the superfluid is formed, slows down the heating or favors the superfluid creation so as to 
weaken the condition (g) effectively by as much as one order of magnitude T{t = Tf orm /10) < 
T c . This amounts to multiplying tl by a factor of 10, giving, under the same conditions as 
above a range of minimal values \a\/r ~ 0.106 — 0.092. 

Finally, we claim that our estimates are lower bounds. In a spherically symmetric 
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T h,k=0 
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harmonic trap the same reasoning as for the homogeneous mixture leads to the follow- 
ing time-dependence for the temperature as a result of background scattering: T/Tp = 
[T /T Ft0 ](N/N )-V«y/l + (3/7r 2 )(7V7o) 2 [l - (AW) 1 / 3 ] 0, with N/N = exp(-t/r L ). 
On the other hand, two- and three-body rates are highest in the high-density middle of the 
trap where the energy per particle is lowest. Thus, lower energy atoms are removed pref- 
erentially fll5| ) giving a temperature that increases more rapidly than in the homogeneous 
system. Secondly, rates for some recombination processes are sensitive to the value of the 
scattering length. A large value of \a\ can imply a shorter lifetime, in which case in fll2|) 



really depends on \a\ [lq|, [17 



In conclusion, we have determined how a metastable degenerate fermion gas heats from 
the hole creation caused by the loss processes. The heating rate is particularly significant at 
ultra-cold temperatures, with a temperature doubling time, r 2 ~ 20[T /T f ] 2 t l . On longer 
time scales t ~ 0.5t£, the fermion system heats up to T > 0.25Tp. For the quest of atom- 
trap fermion superfluidity, this heating competes with the formation of the superfluid and 
we derive a minimal value for the unlike-fermion scattering length to observe superfluid 
formation in a gas mixture of unlike fermion atoms. In contrast to the simple mixture of 
fermion atoms, sympathetic cooling can provide a third (bosonic) atom species that continues 
to absorb the energy released by hole creation. Given the importance of the hole-heating 
mechanism, such continuous cooling may be helpful to reach fermion superfluidity. 

In a question at a Los Alamos colloquium, Albert Pecheck remarked upon the problem of 
Fermi-holes. The fact that this problem did not appear to have been addressed, motivated 
this work. The author also gratefully acknowledges discussions on this subject with Peter 
Milonni, Kyoko Furuya, Randy Hulet, Xinxin Zhao and Michael Di Rosa. 
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Fig. 1 Plot illustrating the loss-induced heating of the fermion mixture. For a homoge- 
neous fermion gas mixture with initial temperature of T = 0.017V, we plot the ratio of 
temperature to fermi-temperature as a function of time for loss processes that are predom- 
inantly background (solid line), two-body (dashed line) and three-body (long-dashed line) 
processes ||13|| . The thin line curve shows T/Tp for the short time temperature approxima- 



tion, T = T Jl + 3t/r 2 in the case of background scattering. 
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Fermi hole heating 




Time (in units of the lifetime x L ) 



